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Abstract. The observed SEDs of a sample of 60 OH/IR stars are fitted using a radiative transfer model of a dusty 
envelope. Among the whole sample, 21 stars have reliable phase-lag distances while the others have less accurate 
distances. L* — P,M — P and M — L, relations have been plotted for these stars. It is found that type E (with 
emission feature at lOp) and type A (with absorption feature at lOp) OH/IR stars have different L* — P and 
M — L* relations while both of them follow a single M — P relation. The type E stars are proven to be located 
in the area without large scale dense interstellar medium while the type A stars are located probably in dense 
interstellar medium. It is argued here that this may indicate the two types of OH/IR stars have different chemical 
composition or zero age main sequence mass and so evolve in different ways. This conclusion has reinforced the 
argument by [Chen et al.(200lj| who reached a similar conclusion from the galactic distribution of about 1000 
OH/IR stars with the IRAS low-resolution spectra (LRS). 
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1. Introduction 


Recently Chen et al.(2001) have published a catalogue of 
OH/IR stars with IRAS LRS spectra. They also argued 
that many type E OH/IR stars which are located at high 
galactic latitude will never evolve into type A OH/IR 
stars . In order to examine the differences between the 
two groups of OH/IR stars (type E and A), it is tried 
in this paper to investigate the luminosity-period relation 
and mass loss rate-period relation of the two groups by 
radiative transfer modeling. The observed data and the 
radiative transfer model used are described in sections 2 
and 3, and a detailed discussion of the derived quantities 
and conclusions can be found in section 4. 


2. Data 

A sample of 60 OH/IR stars with known distances is col¬ 
lected from the literature; 36 of them belong to type E, 18 
belong to type A. For the purpose of comparison, 3 type 
H stars (red continuum), 2 type U stars (unusual spec¬ 
trum) and 1 type F star (featureless) are also included. 
The observational SED (Spectral Energy Distribution) of 
each star consists of three parts: near infrared (NIR) pho¬ 
tometry taken from ground based observations, IRAS far 
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infrared (FIR) photometry and IRAS LRS spectrum. The 
LRS spectra are from the on-line IRAS Data Analysis 
Facility of the University of Calgary. The FIR fluxes are 
obtained from the CD-edition of the published IRAS PSC 
catalogue and only those fluxes with good quality (Q > 2) 
are used. The NIR photometry data of 21 OH/IR stars are 
given by Xiong et al.(1994) and that of the remaining 39 
stars are taken from different papers. The references for 
NIR photometry data are listed in Table 1. The NIR data 
are corrected for interstellar extinction. The formula of 
van Herk(1965)|is used: 


Ay = 0.14csc |6| [1 — exp (—lOdsin | 


( 1 ) 


where b is the galactic latitude, d is the distance from the 
sun. We also use the method from Cardelli et al.(1989) to 


calibrate the data at different wavelengths. 

The periods of pulsation for 46 stars are also searched 
from Chen et al.(2001)| Twinty one sample stars have the 
so called phase-lag distances (this is the most accurate dis¬ 
tance available for OH/IR stars up till now, and hereafter 
these stars are called as the phase-lag sample, and the 
whole 60 stars as whole sample), and these distances are 
taken from Baud and Habing(I983)| Herman et al.(1986)| 


and van Langevelde et al.(1990)|| The distances of the 
other stars are taken from different papers where they are 
estimated by other methods, e.g. kinematic method, us¬ 
ing Loh — B? relation and so on, and so they are not very 
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accurate (error can be larger than 100%). The references 
for distance are listed also in Table 1. For the phase-lag 
sample, the periods, expansion velocities and phase-lag 
distances are shown in Table 2. 


3. Modeling 

3.1. Model description 

The radiative transfer model of a spherical dusty enve¬ 
lope is used to fit the observed SEDs. The model used is 


10 ® 


10 ® 


originally compiled by Szczerba et al.(1997) . This code is 
able to calculate the emerging flux spectrum of a spherical 
shell observed at a specified distance when the luminosity 
(L*) and effective temperature (Teff) of the central star 
and dust species, dust condensation temperature {Tcond), 
dust density distribution {pd ~ r”"), grain size distribu¬ 
tion and the geometrical sizes of the envelope are given. 
The original code is slightly modified so as to make it 
possible to vary only the following five quantities when 
fitting SEDs: Tgff ,Tcond,Tv,n and Ymax, while the other 
quantities are fixed at their typical values: astronomi¬ 
cal silicate grains whose electric permittivity and mag¬ 
netic permeability are taken from [Drain and Lee(1984)| 
and Drain(1985)| grain size distribution of interstellar 
dust as n (a) ~ a~^'^ with the grain size range from 
a,min = 0.005/r to amax = 0.25/ |Mathis et al.(1977) . Here 
Tv is the total optical thickness at 550nm, n is the index 
of negative exponential density distribution law and Ymax 
is the relative geometrical thickness Rmax/Rmin- The ex¬ 
tinction coefficient and the albedo of dust grains are calcu¬ 
lated from Mie theory. In order to fit the observed SEDs, 
2475 models have been calculated by varying Teff from 
1300K to 3800K,Tcond from 800K to 1000K,Ty from 0.1 
to 80,n form 1.4 to 2.6 and Ymax from 10^ to 10^. 

3.2. Model fitting 

The models calculated above are used to fit the ob¬ 
served SEDs simply by weighted linear interpolation. The 
weights are given larger for LRS spectrum and FIR data 
and smaller for the NIR data because the NIR fluxes 
are varied for most sample stars. After having derived 
Teff ,Tv,n,Ymax and the best fitted model SED for each 
star from above, it is easy to determine the luminosity 
(T*) and gas mass loss rate (M) (when the outflow veloc¬ 
ity is known) by assuming the dust to gas mass ratio to 
be 10“^ for all stars. Here it should be pointed out that 
the mass loss rate is calculated at the inner border of the 
circumstellar shell, hence this quantity is theoretically in 
the simultaneous measure with the luminosity of the cen¬ 
ter star. Some derived quantities of the phase-lag sample 
stars are listed also in Table 2. 


4. Analysis on the fitted results 

Firstly the luminosity-period relation (T* —P relation) for 
the phase-lag sample of OH/IR stars is shown in Fig. 1. 
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Fig. 1. The L * —P relation for the phase-lag sample. The 
luminosity shown is obtained by modeling. 

As we can see, for both type E and A stars, the luminosity 
seems to be larger for longer period, but the luminosity 
for the two types seem to be going along parallel paths, 
that is to say, the two T* — P relations do not join with 
each other smoothly from that of the type E to that of 
the type A as we expected. But it is difficult to draw any 
convincing conclusion from Fig. 1, because too few type 
E OH/IR stars are included in this figure. Then the same 
figure for the whole sample is displayed in Fig. 2 where 
many more type E stars are included. This figure is proven 
to confirm the conclusion of Fig. 1. The T* — P relation 
does follow different tracks for the type E and A OH/IR 
stars, although the dispersion in the T* — P relation is 
large for both types. In both Fig. 1 and 2, type H stars 
are located in type A region and seem to follow the T* — P 
relation of type A stars. 

The Ltf—P relation of OH/IR stars and Miras had been 
addressed by many other authors Whitelock et al.(1991)[ 


Le Bertre(1993) Blommaert et al.(I998)| but all of theni 
had considered only a small number of OH/IR stars and 
taken them as a single class. Here we divid the OH/IR 
stars into different subclasses according to their IRAS LRS 
classification and more stars have been included. This give 
us the chance to see some inner structures of the T* — P 
relation for OH/IR stars, as described above. 

The mass loss rate-period relation {M — P relation) 
of the whole sample is plotted in Fig. 3. It is easy to see 
that the whole sample follows a single M — P relation 
and, just as we already know, the mass loss rates increase 
from type E stars with shorter period to type A stars 
with longer period. Similar M — P relation had ever been 
found before for M-type stars by some other authors, eg. 


Whitelock et al.(1994) It should be noted that one star is 
located far from the majority; it is IRAS 18304-0720. This 
star is found to be seriously contaminated by cirrus emis¬ 
sion, and extended structures have been detected around 
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Fig. 2. The L * —P relation for the whole sample. The Fig. 4. The M — L* relation for the OH/IR stars with 


luminosity shown is obtained by modeling. 


known outflow velocity of the whole sample. The mass 
loss rate and luminosity shown are obtained by modeling. 


10 ® 


lO'' 


10 ® 


10 “ 


10 ’ 


10 * 


° S • • 

Cl 


IRAS18304-0720 ^ O 


n 


10^ 


10 ® 

P[days] 


E 

A 

H 


10 '’ 


Fig. 3. The M — P relation for the OH/IR stars with 
known outflow velocity of the whole sample. The mass 
loss rate shown is obtained by modeling. 


it in all four IRAS bands (12, 25, 60,100/r). Hence the ob¬ 
served spectral energy distribution of this star no longer 
represents that of a single shell. 

The M — L* relation is shown in Fig. 4. In this hgure 
the type E stars seem to distribute in an area different 
from that of the stars of other types. Although the lumi¬ 
nosity of type E stars can be very large, the mass loss rate 
seems always to be smaller than that of the other types 
with the same luminosity. 

From the discussion above, the L* — P and M — L* 
relation of type E and A stars seems to follow different 
tracks while the M — P relation of all types in the whole 
sample follow the same track. What do these facts mean 


to the evolution of OH/IR star? According to the general 
evolutionary picture of AGB star, the luminosity and the 
mass loss rate of the AGB stars should increase with time 
and both of them should be the function of the zero age 
main sequence mass (Mzams) and chemical abundance. If 
we assume that all OH/IR stars are evolved from the pro¬ 
genitors of similar M^ams and chemical abundance and 
the type E ones are the precursors of the type A ones as 
we expected, it is reasonable to expect that both the lu¬ 
minosity and the mass loss rate of the type E ones should 
be smaller than that of the type A ones. But in Fig. 4 of 
our paper, many type E stars can reach high luminosity 
while their mass loss rates are always smaller than that of 
the type A stars with the same luminosity. From our Fig. 
1,2,6 and 7, we can also see that some type E stars can 
evolve into luminous objects even if their period is shorter 
than type A ones. This give us the hint that some of the 
two types of OH/IR stars may not evolve from the same 
population of progenitors (with similar Mzams and chem¬ 
ical abundance). If this is the right, we can safely infer 
from the single M — P relation obeyed by the whole sam¬ 
ple that the differences in Mzams and chemical abundance 
of the progenitors have little effect in the M — P relation 
of OH/IR stars. 

Chen et al.(2001) 

have argued that there are probably many type E 
OH/IR stars which would never evolve into type A ones, 
because they found that type A OH/IR stars are mainly 
located on the galactic plane while many type E stars 
are located at high galactic latitude (see their Fig. 5). 
Here a similar figure is given in Fig. 5a which shows a 
similar distribution of different types of OH/IR stars. In 
addition, the more detailed position information of these 
stars are also shown in Fig. 5b and 5c. The distances from 
the stars to the galactic plane and from the projection 
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Fig. 5. a) The space distribution of the whole sample on 
the galactic coordinate, b) Distances of stars to the galac¬ 
tic plane and the galactic center, h: distance to galactic 
plane; r_p: distance from the projections of stars on the 
Galactic plane to the galactic center, c) the projections 
of sample stars on the galactic plane. The galactic center 
is marked by the position of the sun is marked out 
by a dotted line. Numbers beside those circles denote the 


of these stars on the galactic plane to the galactic cen¬ 
ter are plotted in Fig. 5b while the projection of these 
stars on the galactic plane is shown in Fig. 5c. It is indi¬ 
cated from Fig. 5 that the majority of type E sample stars 
are nearby sources while the other types distribute over 
a larger area. There are three galactic arms near the sun: 
the one at several dozen parsecs outside the solar circle is 
the Orion arm; the one at about 2 kpc outside the solar 
circle is the Perseus arm and another at about 2 kpc in¬ 
side the solar circle is the Sagittarius arm. From Fig. 5c, 
it is reasonable to suppose that many type A sample stars 
are located in the three galactic arms. Those stars whose 
distance to the galactic center is smaller than about 5kpc 
must be located inside the dense galactic HI cloud disk. 
And those type E stars are located either in the inter-arms 
area on the galactic plane like the sun or outside of the 
galactic plane. From the space distribution of the sample 
stars, it is not difficult to conclude that, in the whole sam¬ 
ple, the type E and other type stars are probably formed 
in different interstellar environment, hence their chemical 
composition and evolution process can be different. This 
probably is the original cause of the discrepancy between 
the L* — P and M — L* relations of the type E and A 
stars. It should be noted that the distance of many stars 
are not very accurate, so Fig. 5b,c are only valuable for 
reference. 

The conclusions above are mainly based on the quan¬ 
tities derived by modeling. But there exist many uncer¬ 
tainties in the modeling method. The observed data are 
not very certain because of the variation of light, inter¬ 
stellar extinction and emission. The distances to many 
OH/IR stars are not very accurate. There are also many 
uncertainties in the existing radiative transfer modeling. 
For example, the uncertain dust optical properties, the un¬ 
known dust grain physical status and density distribution, 
etc. The spectrum of the center star should be calculated 
in more detail by stellar photosphere models instead of 
using a blackbody spectrum. And many other physical 
processes such as the dust formation process, the dust-gas 
interaction process etc. need to be included in the model. 

Because of these drawbacks of the modeling method, it 
is reasonable to question on the reliability of the derived 
data such as luminosity and mass loss rate. In order to 
clarify the doubts about the derived data, we also calcu¬ 
lated the luminosity by directly integrating the observed 
SEDs. When the spectrum is integrated, it is represented 
by the blackbody spectrum of temperature of 200 K (peak 
at about 25^) in the wavelength range shorter than the 
first NIR observed data point and longer than 100/r. The 
resulted L* — P relations are shown in Fig. 6 for the phase- 
lag sample and in Fig. 7 for the whole sample. In these two 
figures, the conclusion is basically the same. This ensures 
in some degree the reliability of the data derived by mod¬ 
eling and the conclusions derived from them. 

Another problem should also be noted: some of the 
sample stars, even some type E ones, have their luminos¬ 
ity very large (about IO^Lq); this is unusual for OH/IR 
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Fig. 6. The L* — P relation for the phase-lag sample. The 
luminosity shown is obtained by directly integrating the 
observed SEDs. 
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Fig. 7. The L* — P relation for the whole sample. The 
luminosity shown is obtained by directly integrating the 
observed SEDs. 


stars. In order to verify whether any calculation errors 
have occurred, the luminosity data of some sample stars 


are compared with that from van Langevelde et al.(1990) 


who had estimated luminosity for some OH/IR stars with 
phase-lag distances in another way. The luminosity data 
derived in this paper are found to be comparable with 
theirs. And there are also some stars in their paper reach a 
high luminosity near to IO^Pq. It is reasonable to suspect 
that these luminous stars may actually be supergiants. 
But further inspection is needed for this subject. 


used in our work the online SIMBAD database and the online 
IRAS Data Analysis Facility at the University of Calgary. This 
work is supported by grants from National Natural Science 
Foundation of China and Chinese Academy of Sciences. 
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Table 1. References for NIR and distance data. 


IRASjiame NIR dist. 

IRASjiame NIR dist. 

IRAS name 

NIR dist. 

01037+1219 

21 

15 

14297-6010 

15 

13 

23412-1533 

5 

13 

01304+6211 

21 

15 

15193+3132 

22 

14 

23425+4338 

22 

16 

03293+6010 

18 

3 

16474-4418 

15 

13 

15255+1944 

21 

1 

03336-7636 

5 

16 

16494-4327 

15 

13 

18182-1504 

7 

9 

03507+1115 

5 

9 

17073-4225 

15 

13 

18257-1000 

2,4,7 

19 

04140-8158 

5 

16 

17073-3955 

15 

13 

18257-1052 

2,11,17 

19 

05098-6422 

5 

16 

17271-3425 

15 

13 

18276-1431 

15 

9,1 

05528+2010 

21 

1 

17317-3331 

15 

13 

18348-0526 

21 

19 

06500+0829 

12 

16 

17360-3012 

15 

13 

18431-0403 

8,15 

19 

07209-2540 

5 

13 

17540-1919 

12 

14 

18432-0149 

8,15,17 

9 

07585-1242 

5 

16 

18135-1456 

15 

9 

18445-0238 

4,11,15 

9 

09425+3444 

22 

16 

18437-0643 

15 

9 

18460-0254 

2,4,11,15 

9 

09429-2148 

5 

3 

19017+0608 2,18 

3 

18488-0107 

2,8,10,11 

19 

10189-3432 

5 

16 

19059-2219 

5 

13 

18498-0017 4,11,12,17 

19 

10287-5733 

15 

13 

20000+4954 

18 

9 

18549+0208 

8 

19 

10580-1803 

5 

16 

20077-0625 

5 

13 

18560+0638 

21 

19 

12310-6233 

15 

13 

20259-4035 

20 

14 

19039+0809 

21 

9 

13157-6421 

15 

13 

21206-4054 

20 

14 

19192+0922 

6 

19 

14086-2839 

5 

14 

22516+0838 

5 

13 

19550-0201 

5 

1 

14247+0454 

5 

9 

23041+1016 

5 

14 

22177+5936 

21 

9 


Refer ence list: (l)|Baud and Habing(198'^ (2) [Baud et al.(1985)] ( 3) [Bowers et al.(1983)| (4 ) [Evans and Bechwith(1984) 
(5) [Fouque et al.(1992)[ (6) [Fix and Mutel(1984)[ (7) [Garcia-Lario et al.( 1997)1 (8) [Gehrz et al.(198^ 
(9) [Herman et al.(1986'^ (1 0) [Kwok et al.( 1987) [ (11) [Jones et al.(1983)[ (12) "* [Lepine et al.(1995)[ (13) 
[Le Sidaner and Le Bertre(1996)| (14) [Nguyen-Q-Rie et al.(1979j| (15) [Nyman et al.(1993)[ (16) [Ortiz and Maciel( 1996) [ 
(17) [Schultz et al.(1976)[ (18) Sun and Zhang(1998)[ (19) [van Langevelde et al.(1990)[ (20) [Whitelock et al.(1994)[ (21) 
[Xiong et al.(1994)[ (22) query from SIMBAD database 


Table 2. Some parameters of the stars with phase-lag distances. 


IRASjiame 

period 14 

[day] [km/s] 

D 

[kpc] 

L* M LRS type 

[lO^L©] [Mg/yr] 

01037+1219.... 

.. 660.0 

18.2 

0.74 

24.56 

0.96 

E 

05528+2010.... 

.. 368.3 

2.7 

1.20 

76.86 

0.11 

E 

15255+1944.... 

.. 425.1 

7.5 

3.10 

105.71 

0.58 

E 

19039+0809.... 

.. 284.2 

6.6 

0.29 

3.46 

0.02 

E 

18431-0403.... 

.. 717.0 

17.5 

1.45 

8.36 

0.24 

E 

19550-0201.... 

.. 394.8 

5.4 

0.85 

27.20 

0.06 

E 

01304+6211.... 

.. 1995.0 

11.0 

2.90 

50.91 

6.86 

A 

18348-0526.... 

.. 1566.0 

14.0 

1.44 

18.09 

4.50 

A 

18560+0638.... 

.. 1424.0 

16.1 

2.04 

22.31 

4.93 

A 

22177+5936.... 

.. 1215.0 

14.5 

2.38 

14.54 

5.58 

A 

18182-1504.... 

.. 1011.7 

20.5 

0.55 

0.84 

0.78 

A 

18257-1052.... 

.. 1129.5 

18.2 

5.06 

13.09 

2.04 

A 

18257-1000.... 

.. 1975.2 

19.0 

8.51 

83.31 

18.99 

A 

18276-1431.... 

.. 890.0 

12.0 

2.08 

3.09 

0.60 

H 

18432-0149.... 

.. 1038.6 

17.5 

7.70 

40.99 

8.26 

A 

18445-0238.... 

.. 853.0 

17.0 

1.09 

1.46 

0.58 

A 

18460-0254.... 

.. 2064.4 

19.5 

1.77 

9.67 

3.59 

A 

18488-0107.... 

.. 1540.2 

20.6 

11.87 

81.31 

16.32 

A 

18498-0017.... 

.. 1536.3 

15.9 

5.02 

22.58 

1.58 

H 

18549+0208.... 

.. 840.0 

13.6 

4.21 

6.30 

0.91 

A 

19192+0922.... 

.. 552.0 

16.1 

1.13 

13.82 

0.52 

A 




















































